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Abstract 


DRDC  Atlantic  has  developed  a  coherent  transmission-loss  model  that  simulates 
the  propagation  of  acoustic  pulses  through  an  ocean  environment.  Using  a  set  of 
input  eigenrays  and  an  input  waveform,  the  WATTCH  (Waveform  Transmission 
Through  a  Channel)  model  can  simulate  the  signal  received  as  a  result  of  trans¬ 
mitting  the  waveform  through  the  ocean  environment.  Multiple  output  time-series 
channels  represent  receivers  at  given  ranges  and  depths.  Assuming  the  required 
eigenrays  can  be  generated,  WATTCH  can  simulate  the  effects  of  complex  range- 
dependant  environments. 

This  document  presents  the  mathematical  formulation  and  set  of  examples  used  to 
develop  and  verify  the  model.  In  addition,  a  comparison  is  made  between  two  pulse 
propagation  techniques.  The  first  technique  uses  WATTCH  to  simulate  the  arrival 
of  a  transmitted  waveform.  The  second  technique  uses  WATTCH  to  simulate  the 
arrival  of  a  band-limited  impulse  waveform,  and  convolves  the  results  with  the  de¬ 
sired  transmitted  waveform.  The  comparison  shows  the  techniques  yield  equivalent 
results  and  therefore  multiple  waveforms  may  be  simulated  using  the  second  tech¬ 
nique,  but  only  running  the  WATTCH  model  once. 


Resume 


RDDC  Atlantique  a  elabore  un  modele  coherent  d’affaiblissement  de  transmis¬ 
sion  qui  simule  la  propagation  d’ impulsions  acoustiques  dans  un  environnement 
oceanique.  En  utilisant  en  entree  un  ensemble  de  rayons  propres  et  une  forme 
d’onde,  le  modele  WATTCH  (Waveform  Transmission  Through  a  Channel,  trans¬ 
mission  d’une  forme  d’onde  dans  un  canal)  peut  generer  un  ensemble  de  serer 
series  chronologiques.  Chaque  serie  chronologique  represente  les  signaux  prevus 
correspondant  une  distance  et  une  profondeur  donnees.  Si  les  rayons  propres  re- 
quis  peuvent  Itre  generes,  WATTCH  peut  simuler  les  effets  d’environnements  com¬ 
plexes. 

Le  present  document  montre  la  formulation  mathematique  et  un  ensemble  d’exemples 
utilises  pour  elaborer  et  verifier  le  modele.  En  outre,  il  montre  une  comparaison 
entre  deux  techniques  de  propagation  d’impulsions.  La  premiere  technique  utilise 
WATTCH  pour  simuler  la  reception  d’une  forme  d’onde  transmise.  La  deuxieme 
technique  utilise  WATTCH  pour  simuler  la  reception  d’une  forme  d’ impulsion 
bande  limitee,  puis  elle  convolutionne  les  resultats  avec  la  forme  d’onde  transmise 
voulue.  La  comparaison  montre  que  les  techniques  donnent  des  resultats  equivalents, 
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de  sorte  que  plusieurs  formes  d’onde  peuvent  itre  simulees  l’aide  de  la  deuxieme 
technique,  mais  en  necessitant  une  seule  execution  du  modele  WATTCH. 

Ceci  est  le  resume  en  fran<jais. 
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Executive  summary 

Background 

Underwater  acoustic  communication  requires  the  robust  transfer  of  information 
through  the  water  column.  Communication  systems  must  convert  a  sequence  of 
digital  signals  to  a  sequence  of  acoustic  symbols  that  are  transmitted  through  the 
underwater  environment.  A  receiving  system  must  convert  the  acoustic  signals  in 
order  to  process  and  reconstruct  the  original  digital  signals. 

The  underwater  acoustic  environment  impacts  on  the  received  signals.  Multiple 
arrivals  and  signal  distortion  may  cause  degradation  in  the  arrived  signals.  This 
degradation  may  result  in  a  loss  of  ability  to  reconstruct  the  original  information. 


Results 

The  WATTCH  model  predictions  are  compared  with  other  methods  of  estimating 
the  impact  of  the  acoustic  environment  on  a  communication  system.  The  WATTCH 
model  is  then  used  to  produce  simulated  received  signals  based  on  two  approaches. 
The  first  approach  uses  WATTCH  to  simulate  the  received  time  series  based  on 
transmitting  a  given  waveform.  The  second  approach  uses  WATTCH  to  simulate 
the  received  time  series  based  on  transmitting  a  band-limited  impulse  waveform. 
The  output  is  convolved  with  the  original  waveform  to  yield  a  received  time  series 
equivalent  to  the  first  approach. 


Significance 

The  WATTCH  model  has  been  delivered  for  use  in  signal  processing  studies.  The 
impact  of  the  environment  on  underwater  acoustic  communication  systems  is  of 
particular  interest,  but  the  model  can  be  used  to  investigate  coherent  pulse  propaga¬ 
tion  for  Anti-Submarine  Warfare  active  sonar  studies. 


Future  work 

The  WATTCH  model  is  to  be  used  to  study  underwater  acoustic  communication. 
Extensions  are  likely  to  be  undertaken  in  order  to  interpret  measurements. 


James  A.  Theriault  and  Sean  Pecknold;  2006;  Impulse  Propagation 
using  WATTCH;  DRDC  Atlantic  ECR  2004-248;  Defence  R  &  D 
Canada  -  Atlantic. 
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Sommaire 


Introduction 

Les  communications  acoustiques  sous-marines  necessitent  un  transfert  fiable  d’ in¬ 
formation  dans  la  colonne  d’eau.  Les  systemes  de  communication  doivent  convertir 
une  sequence  de  signaux  numeriques  en  une  sequence  de  symboles  acoustiques  qui 
sont  transmis  dans  l’environnement  sous-marin.  Un  systeme  de  reception  doit  effec- 
tuer  la  transduction  des  signaux  acoustiques  pour  traiter  et  reconstruire  les  signaux 
numeriques  initiaux. 

L’environnement  acoustique  sous-marin  a  une  incidence  sur  les  signaux  reus.  Les 
receptions  multiples  et  la  distorsion  des  signaux  peuvent  causer  une  degradation  des 
signaux  reus,  laquelle  peut  reduire  la  capacite  de  reconstruction  de  l’information 
initiale. 

Resultats 

Les  previsions  du  modele  WATTCH  sont  comparees  aux  previsions  obtenues  avec 
d’autres  methodes  d’estimation  de  l’incidence  de  l’environnement  acoustique  sur 
un  systeme  de  communication.  Le  modele  WATTCH  est  alors  utilise  pour  produire 
deux  signaux  reus,  le  premier  base  sur  l’hypothese  d’une  forme  d’onde  donnee,  le 
deuxieme  utilisant  une  forme  d’ impulsion  bande  limitee.  En  convolutionnant  les 
resultats  de  la  prevision  obtenue  pour  une  forme  d’ impulsion  avec  la  forme  d’onde 
donnee,  on  montre  que  le  resultat  est  equivalent  la  transmission  simulee  de  la  forme 
d’onde  initiale. 

Portee 

Le  modele  WATTCH  a  ete  livre  en  vue  de  son  utilisation  dans  les  etudes  de  trai- 
tement  des  signaux.  L’incidence  de  l’environnement  sur  les  systemes  de  commu¬ 
nications  acoustiques  sous-marines  presente  un  interlt  particulier,  mais  le  modele 
peut  Itre  utilise  pour  etudier  la  propagation  coherente  d’ impulsions  aux  fins  de  son 
application  au  sonar  actif  de  guerre  anti-sous-marine. 


Recherches  futures 

Le  modele  WATTCH  est  destine  a  etre  utilise  pour  T etude  des  communications 
acoustiques  sous-marines  et  il  sera  elargi  seulement  a  la  suite  d’une  interpretation 
plus  approfondie  des  mesures  effectuees  dans  l’essai  en  mer. 

Ceci  est  le  sommaire  en  frangais. 


IV 
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1  Introduction 


Underwater  acoustic  communication  requires  the  robust  transfer  of  information 
through  the  environment.  Systems  for  communication  must  convert  a  sequence 
of  digital  signals  to  a  sequence  of  acoustic  symbols  that  are  transmitted.  A  receiv¬ 
ing  system  must  convert  the  acoustic  signals  in  order  to  process  and  reconstruct  the 
original  digital  signals. 

The  underwater  acoustic  environment  impacts  on  the  received  signals.  Multiple 
arrivals  and  signal  distortion  may  cause  degradation  in  the  arrived  signals.  This 
degradation  may  result  in  a  loss  of  ability  to  reconstruct  the  original  information. 

With  support  from  the  US  Office  of  Naval  Sonar  (ONR),  DRDC  Atlantic  has  devel¬ 
oped  a  coherent  pulse-propagation  model.  Though  full-wave-equation  solutions, 
such  as  OASES  [1],  exist,  they  are  computationally  intensive  and  restricted  in  the 
environments  they  can  model.  Westwood  [2]  claimed  that  a  ray-theory  based  model 
could  contain  all  of  the  relevant  physics  for  the  problem  of  interest,  and  be  signifi¬ 
cantly  less  computationally  demanding. 

Starting  with  a  set  of  input  eigenrays  and  an  input  waveform,  the  WATTCH  (Wave¬ 
form  Transmission  Through  a  Channel)  model  can  simulate  the  expected  time  series 
received  at  a  set  of  hydrophones.  Each  time  series  represents  the  expected  signals 
corresponding  to  a  given  range  and  depth.  Assuming  the  required  eigenrays  can  be 
generated,  WATTCH  can  simulate  the  effects  of  complex  range -dependant  envi¬ 
ronments. 

This  document  presents  the  mathematical  formulation  and  set  of  examples  used  to 
develop  and  verify  the  model.  For  the  purposes  of  this  document,  the  eigenrays 
are  generated  by  the  US  Generic  Sonar  Model  (GSM)  [4].  However,  in  principle, 
eigenrays  generated  from  other  models  such  as  CASS/GRAB  [5]  or  Bellhop  [6] 
could  also  be  used.  The  results,  based  on  GSM,  are  compared  with  time  series 
predicted  by  Chapman  et  al.  [7]  using  the  benchmark  OASES  model.  A  separate 
document  [8]  describes  the  software.  Using  the  benchmark  model  (OASES)  [1]  at 
a  2  km  range  on  a  dual  2  GHz  Xeon  computer  (with  2  GB  of  memory)  running 
Linux,  10.44  hours  were  required.  The  same  test,  using  GSM/WATTCH  required 
9  s.  These  times  and  the  agreement  between  the  benchmark  and  WATTCH  model 
support  Westwood’s  conclusion  [2]. 

In  addition,  a  comparison  is  made  between  two  pulse  propagation  techniques.  The 
first  technique  uses  WATTCH  to  simulate  the  arrival  of  a  transmitted  waveform. 
The  second  technique  uses  WATTCH  to  simulate  the  arrival  of  a  band-limited  im¬ 
pulse  waveform,  then  convolve  the  time  series  with  the  desired  transmitted  wave¬ 
form. 
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Using  a  waveform  and  environment  presented  by  Chapman  et  al.  [7],  together  with 
a  500  Hz  to  6500  Hz  band-limited  impulse  waveform,  the  comparison  shows  that 
the  techniques  yield  equivalent  results.  Therefore  multiple  waveforms  may  be  sim¬ 
ulated  using  the  second  technique  while  only  running  the  WATTCH  model  once. 
The  advantage  with  the  second  technique  is  that  the  time  required  for  execution  of 
the  WATTCH  model  depends  on  the  sampling  frequency  and  length  of  the  modeled 
waveform.  The  second  technique  involves  propagating  a  very  short  duration  wave¬ 
form,  and  will  generally  require  significantly  less  time  for  any  reasonable  wave¬ 
form  duration.  In  the  example  given,  the  WATTCH  model  required  nine  seconds 
to  propagate  the  band-limited  impulse  waveform  to  a  distance  of  2  km,  while  a 
0.2  second  long  LFM  pulse  required  33  seconds  of  computation.  The  convolution 
required  approximately  15  s.  Therefore,  by  using  the  convolution  approach,  ap¬ 
proximately  9  s  of  computational  load  are  saved  for  each  receiver  distance  for  this 
waveform,  with  more  time  saved  for  any  additional  waveform  types  required.  This 
computational  savings  is  increased  with  increasing  waveform  length  and  sampling 
rate.  Additionally,  in  the  case  where  the  pulse  to  be  propagated  is  either  very  long 
or  of  very  large  bandwidth,  the  WATTCH  model  can  become  memory-limited.  This 
problem  is  eliminated  by  propagating  the  impulse  function. 
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2  The  WATTCH  Model 


The  WATTCH  model  simulates  the  effect  of  transmitting  a  waveform  through  an 
underwater  acoustic  environment.  Figure  1  shows  the  basic  operations  required  to 
produce  a  simulated  time  series.  Starting  with  a  description  of  the  acoustic  environ¬ 
ment,  geometry,  and  frequency  band,  a  set  of  eigenrays  is  generated.  The  eigenrays 
represent  all  of  the  significantly-contributing  acoustic  paths  connecting  the  source 
to  a  receiver.  The  source  location  is  assumed  to  be  defined  by  a  depth,  zq-  at  range 
zero.  For  a  given  receiver  depth,  z,  the  range,  r  can  vary  by  equal  increments  be¬ 
tween  a  minimum  and  maximum.  The  frequency  range  is  assumed  to  cover  the 
spectrum  of  the  intended  waveforms. 

The  eigenrays  can  be  computed  by  various  acoustic  models  such  as  the  US  Generic 
Sonar  Model  (GSM)  [4],  CASS/GRAB  [5],  or  Bellhop  [6].  GSM  is  used  as  the 
baseline  model.  However,  it  was  recognized  at  the  outset  that  the  enhanced  ability 
for  either  GRAB  or  BellHop  to  model  an  environment  that  changes  with  range  may 
be  required.  Investigations  of  using  both  GRAB  and  BellHop  have  been  carried  out 
as  separate  activities. 

Using  the  set  of  eigenrays  and  the  waveform  (specified  in  DRDC  .dat32  format  [9]), 
WATTCH  computes  the  predicted  time  series.  The  process  starts  with  computing 
the  Fast  Fourier  Transform  (FFT)  of  the  waveform.  This  yields  the  frequency- 
dependent  amplitude  and  phase  representation  of  the  waveform.  Conceptually, 
each  of  the  frequency  components  then  propagates  as  described  by  the  frequency- 
dependant  eigenrays.  The  propagated  components  are  then  reassembled  to  yield 
the  predicted  waveforms  at  the  defined  receivers. 

The  mathematical  formulation  starts  with  the  real  input  waveform,  s(t),  which  can 
be  represented  using  a  discrete  Fourier  series, 

oo 

s(t)  =  K  COs(27 lfnt  +  <| )„)  (1) 

72=0 

where  s  is  defined  for  time  t  between  0  and  To.  the  pulse  length.  The  amplitude,  bn, 
and  phase,  §n  are  the  real  discrete-Fourier  Transform  coefficients  at  the  frequencies 
fn.  The  previous  equation  can  be  approximated 

N 

(2) 

22=0 


fn  is  defined  to  be  n/  (NT )  where  N  and  T  are  the  number  of  points  in  the  FFT  and 
the  sampling  period,  respectively,  and  sn(t)  —  cos(2%f„t  +  §n) . 
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Figure  1:  WATTCH  model  concept 

Consider  propagating  this  waveform  through  the  underwater  acoustic  medium.  Mul¬ 
tiple  paths  connect  the  source  and  receiver  (see  Figure  2).  For  a  given  geometry 
(r.z,Zo  specified),  environment,  and  eigenray  tolerance  (defining  if  an  eigenray  is 
making  a  significant  contribution),  there  are  M  eigenrays.  Each  eigenray  connect¬ 
ing  the  source  and  a  receiver  defines  the  parameters  describing  acoustic  propagation 
conditions.  Each  eigenray  is  parameterized  as  a  function  of  frequency,  /,  by  the 
launch  angle,  arrival  angle,  amplitude  am,  phase  shift  0m,  and  time  of  flight  tm.  For 
WATTCH  purposes,  the  launch  and  arrival  angles  are  not  required.  Though  both 
the  time  of  flight,  tm.  and  the  phase  shift,  0,n  can  be  assumed  to  be  independent  of 
frequency,  only  tm  will  be  assumed  to  be  independent  of  frequency. 

Considering  the  mth  eigenray,  let  amM(r.z,Zo)  and  Qm,n(r-Z-Zo)  represent  the  ampli¬ 
tude  and  phase  shift  at  the  ffj1  frequency.  Hence,  each  of  the  sinusoidal  components 
sn  needs  to  be  shifted  by  Qm,n(r,z,zo)  resulting  in  the  received  time  series  from  the 
pulse  component,  sn(r,z,zo),  along  the  mth  eigenray  to  be  given  by 


Mm,n  (t  jPZ,  Zo) 


am,n(p  Z.- Zo)bn  X 

(cos(27t/„r  +  ()„  Qm  n(r. z.  z.q)) :  if  t  G  \tm .  tm  T  To] 
0  Othem’ise. 


Hence,  the  simulated  time  series,  u,  that  includes  contributions  from  all  of  the  eigen- 
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Figure  2:  Example  multi-path  ray  diagram 

rays  and  frequency  components  is  given  by 

M(r,z,zo)  N 

u(t,r,z,zo)  =  ^2  y^MOT;n(i,r,z,zo)  (4) 

m=  1  /;= 0 


The  WATTCH  model  uses  this  formulation  to  simulate  the  time  series  at  a  given 
receiver.  After  choosing  source  and  receiver  depths,  and  specifying  a  set  of  ranges, 
WATTCH  will  produce  a  time  series  associated  with  each  range.  The  current  ver¬ 
sion  of  WATTCH  adjusts  the  time  scale  by  subtracting  the  minimum  tm(r,z,zo ) 
for  a  given  range  from  the  associated  time  series.  By  shifting  the  time-series,  in 
time,  the  arrival  information  is  preserved  while  requiring  significantly  less  com¬ 
puter  memory  and  disk  space.  The  software  is  described  by  Calnan  [8]. 
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3  WATTCH  Verification 


3.1  Constant  Sound  Speed  Environment 

Section  2  describes  the  WATTCH  mathematical  formulation.  Calnan  [8]  imple¬ 
mented  the  described  model  in  the  IDL  language.  Chapman  et  al  [7]  presented 
the  development  and  comparison  of  a  set  of  benchmark  models.  Starting  with 
OASES  [1]  as  a  benchmark,  Reference  [7]  compared  it  with  the  Mathematica-based 
KosmicRay  model.  Being  a  full- wave  model,  OASES  is  a  theoretically  exact,  ana¬ 
lytic  solution  of  the  wave  equation,  and  is  thus  taken  as  the  benchmark  “truth.” 

Figure  3  shows  the  results  presented  as  Figure  15  in  Reference  [7].  The  two  chan¬ 
nels  shown  in  the  figure  represent  the  simulated  arrivals  of  a  500  to  6500  Hz  band- 
limited  function  at  a  range  of  2  km.  The  pulse  length  is  0.01  s.  The  broadband 
waveform  is  represented  by 


1  ffl 

s(t )  =  - —  /  cos(2;tfk)  JP  (5) 

J2  —  JI  J fx 

where  f\  and  fa  represent  the  low  and  high  frequency  covered  by  the  bandwidth.The 
environment  has  a  120  m  depth,  with  the  source  and  receiver  at  100  m.  The  sound 
speed  is  1450  m/s.  A  Raleigh  phase  shift  [3]  and  reflection  coefficient  is  assumed 
for  the  seabed  reflection.  The  bottom  density  is  1.9  g/cm3  and  the  bottom  sound 
speed  is  1650  m/s.  The  two  models  show  reasonably  good  agreement  (for  example, 
see  the  circled  sections  considered  in  detail  in  Reference  citechapman). 

This  environment  was  modeled  by  GSM  and  used  for  input  in  the  WATTCH 
model.  Alternatively,  KosmicRay  is  used  to  generate  the  eigenray  input.  Figure  4 
shows  the  comparison  between  the  WATTCH  model  with  the  GSM  and  Kosmi¬ 
cRay  eigenray s.  The  results  based  on  the  KosmicRay  input  agree  with  those  gen¬ 
erated  in  Reference  [7],  but  the  results  based  on  the  GSM  input  show  noticeable 
discrepancy.  The  GSM-based  eigenrays  have  different  travel  times  than  those  gen¬ 
erated  by  KosmicRay,  which  result  in  some  of  the  paths  adding  out  of  phase  rather 
than  in  phase. 

This  may  be  a  limitation  of  GSM’s  capabilities  in  calculating  eigenrays  for  iso¬ 
velocity  profiles.  However,  it  is  also  the  case  that  even  slight  changes  in  environ¬ 
ment  may  yield  significant  changes  in  eigenray  travel  times,  phases  and  amplitudes. 
If  the  waveform  is  sensitive  to  these  changes,  the  output  time  series  may  be  altered 
significantly  as  well.  In  the  isovelocity  case,  the  difference  in  eigenray  arrival  times 
calculated  by  GSM  is  sufficient  to  cause  these  changes. 
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Figure  3:  Comparison  of  time  series  generated  by  OASES  and  KosmicRay  for  a 
band-limited  impulse  function  at  2  km  distance  in  a  shallow  isovelocity  ocean 
environment. 

Figure  5  shows  the  outputs  assuming  the  same  waveform  and  environment  used 
above,  with  the  exception  that  in  one  case,  the  surface  sound  speed  has  been  changed 
to  1453.25  m/s  (negative  linear  sound  speed  gradient).  In  the  second  case  the  bot¬ 
tom  sound  speed  has  been  changed  to  1453.25  m/s  (positive  linear  sound  speed 
gradient). 

In  these  cases,  only  the  first  part  of  the  first  group  of  arrivals  shows  a  noticeable 
difference  between  the  two  examples,  and  both  show  a  greater  degree  of  similarity 
towards  the  isovelocity  output  of  KosmicRay  than  to  the  isovelocity  output  of  GSM 
for  most  of  the  arrival  groups.  This  is  consistent  with  previous  experience  with 
GSM.  Isovelocity  profiles  often  cause  difficulty  for  GSM. 


3.2  Bilinear  Sound  Speed  Profile  Environment 

One  of  the  advantages  of  using  GSM  rather  than  KosmicRay  as  the  source  of 
eigenrays  is  that  it  is  not  restricted  to  isovelocity  environments.  An  example  of  this 
is  a  bilinear  profile,  with  a  negative  sound  speed  gradient  near  the  surface  which 
changes  to  a  positive  gradient  for  deeper  depths. 

In  this  example,  the  pulse  used  is  the  same  as  in  the  isovelocity  case,  a  0.01  s 
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Reduced  time  t-r/1 .454  (s) 

Figure  4:  Comparison  of  time  series  generated  by  WATTCH  based  on  eigenrays 
from  GSM  and  KosmicRay  for  band-limited  impulse  function  at  2  km  distance  in  a 
shallow  isovelocity  ocean  environment. 

band-limited  (500  -  6500  Hz)  impulse  function.  The  water  depth  is  120  m,  with 
source  and  receiver  depths  of  100  m.  The  sound  speed  is  bilinear,  with  surface 
speed  of  1480  m/s.  This  changes  linearly  to  1442  m/s  at  60  m  then  increases  lin¬ 
early  to  1447  m/s  at  the  seabed.  An  isovelocity  model  such  as  KosmicRay  cannot 
address  this  problem.  The  bottom  density  is  1.9  g/cm3  and  the  bottom  sound  speed 
is  1650  m/s.  A  Rayleigh  phase  shift  [3]  and  bottom  reflection  coefficient  is  used. 

Figure  6  shows  the  OASES  and  WATTCH  results.  The  two  models  show  generally 
good  agreement  in  this  case. 
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Figure  5:  Comparison  of  time  series  generated  by  WATTCH  based  on  eigen  rays 
from  GSM  for  band-limited  impulse  function  at  a  2  km  distance  in  a  shallow  ocean 
environment,  with  slightly  differing  linear-gradient  sound  speed  profiles. 
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Figure  6:  Comparison  of  time  series  generated  by  WATTCH  and  OASES  for  for  a 
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ocean  environment. 
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4  Broadband  pulse  convolution 


Section  3  verified  the  operation  of  the  WATTCH  model.  It  is  capable  of  producing 
a  simulated  time  series  that  can  be  used  in  studying  the  coherence  limits  of  an 
environment.  In  order  to  save  computational  resources  and  to  reduce  the  effort  for 
the  modeler,  an  alternative  to  computing  the  time  series  from  a  given  waveform  is 
considered. 

Computing  the  received  time  series  based  on  a  given  waveform  may  be  accom¬ 
plished  by  convolving  a  transfer  function  with  the  original  waveform.  The  transfer 
function  may  be  generated  by  simulating  the  propagation  of  a  band- limited  impulse 
function.  Using  this  approach  requires  only  one  execution  of  the  WATTCH  model 
for  a  given  environment  and  set  of  geometries.  The  impact  of  the  environment  on 
a  set  of  waveforms  may  be  studied  without  repeatedly  executing  the  model.  This 
results  in  significant  time  savings. 

This  section  shows  the  equivalence  of  the  two  approaches.  Care  must  be  taken  when 
calculating  the  final  convolution,  as  some  software  packages  may  not  compute  a 
convolution  as  expected. 

Figures  7  shows  a  comparison  between  the  direct  computation  through  WATTCH 
of  a  propagated  LFM  and  the  same  LFM  convolved  with  the  WATTCH-computed 
band-limited  impulse  function  described  above.  The  environment  is  the  bilinear 
shallow  ocean  described  in  Section  3.2.  The  Hanning-windowed  LFM  has  duration 
0.2  s,  with  a  centre  frequency  of  3000  Hz  and  a  bandwidth  of  400  Hz. 

Figure  8  shows  the  power  spectra  for  the  waveform  at  2  km  using  both  approaches. 
As  expected,  the  power  spectrum  for  the  time  series  generated  by  convolving  the 
LFM  with  the  band-limited  impulse  function  shows  significantly  less  power  in  the 
frequency  bands  outside  of  the  500  to  6500  Hz  band  of  the  band-limited  impulse 
function  and  shows  good  agreement  within  the  bandwidth.  Hence,  the  convolution 
approach  works  well  for  waveforms  having  the  majority  of  its  energy  in  the  same 
band  as  the  “impulse  waveform.”  Conversely,  for  waveforms  with  significant  en¬ 
ergy  outside  of  the  “impulse  waveform,”  the  approach  would  introduce  significant 
biases. 
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Figure  7:  Time  series  generated  by  WATTCH  for  an  LFM  (shallow  bilinear-profile 
ocean  environment)  (upper  trace)  and  by  convolving  the  LFM  with  WATTCH 
output  (lower  trace)  based  on  a  band-limited  impulse  waveform  at  different  ranges. 


Figure  8:  Power  spectrum  comparison  for  the  2  km  outputs  for  Figure  7. 
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5  Summary 


Underwater  acoustic  communication  requires  the  robust  transfer  of  information 
through  the  water  column.  Communication  systems  must  convert  a  sequence  of 
digital  signals  to  a  sequence  of  acoustic  symbols  that  are  transmitted  through  the 
underwater  environment.  A  receiving  system  must  convert  the  acoustic  signals  in 
order  to  process  and  reconstruct  the  original  digital  signals. 

The  underwater  acoustic  environment  impacts  on  the  received  signals.  Multiple 
arrivals  and  signal  distortion  may  cause  degradation  in  the  arrived  signals.  This 
degradation  may  result  in  a  loss  of  ability  to  reconstruct  the  original  information. 

The  WATTCH  model  predicts  the  received  time  series  based  on  eigenray  input, 
geometry,  and  waveforms.  This  document  presents  the  formulation,  verification, 
and  an  alternative  approach  to  directly  computing  the  resulting  time  series. 
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